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Abstract: We provide a holographic evaluation of novel contributions to the drag force 
acting on a heavy quark moving through strongly interacting plasma. The new contributions 
are chiral in the sense that they act in opposite directions in plasmas containing an excess 
of left- or right-handed quarks. The new contributions are proportional to the coefficient of 
the axial anomaly, and in this sense also are chiral. These new contributions to the drag 
force act either parallel to or antiparallel to an external magnetic field or to the vorticity of 
the fluid plasma. In all these respects, these contributions to the drag force felt by a heavy 
quark are analogous to the chiral magnetic effect (CME) on light quarks. However, the new 
contribution to the drag force is independent of the electric charge of the heavy quark and 
is the same for heavy quarks and antiquarks, meaning that these novel effects do not in fact 
contribute to the CME current. We show that although the chiral drag force can be non¬ 
vanishing for heavy quarks that are at rest in the local fluid rest frame, it does vanish for 
heavy quarks that are at rest in a suitably chosen frame. In this frame, the heavy quark at 
rest sees counterpropagating momentum and charge currents, both proportional to the axial 
anomaly coefficient, but feels no drag force. This provides strong concrete evidence for the 
absence of dissipation in chiral transport, something that has been predicted previously via 
consideration of symmetries. Along the way to our principal results, we provide a general 
calculation of the corrections to the drag force due to the presence of gradients in the flowing 
fluid in the presence of a nonzero chemical potential. We close with a consequence of our 
result that is at least in principle observable in heavy ion collisions, namely an anticorrelation 
between the direction of the CME current for light quarks in a given event and the direction 
of the kick given to the momentum of all the heavy quarks and antiquarks in that event. 



Contents 


1 Introduction and Summary 1 

2 Holographic fluid 4 

3 Drag Force 7 

3.1 Fluid at rest 7 

3.2 Corrections to the Drag Force Due to Fluid Gradients 10 

4 Anomalous Contributions 12 

4.1 Chiral Vortical Drag Force 13 

4.2 Chiral Magnetic Drag Force 14 

4.3 Contributions to the Drag Force arising from the Gravitational Anomaly 15 


5 The Dissipationless Character of the Chiral Magnetic and Vortical Effects 16 

6 Outlook and Potential Phenomenological Consequences 20 


1 Introduction and Summary 

The analysis of how a heavy quark moving through the strongly coupled liquid quark-gluon 
plasma produced in ultrarelativistic heavy ion collisions loses energy is motivated by heavy 
ion collision experiments, in which these heavy quarks are used as probes of the plasma and, 
for those that lose enough energy, as tracers that follow its flow. If one assumes that the 
interactions between the heavy quark and the plasma are weak, then perturbative methods 
that were first introduced in Ref. [1] can be employed to analyze heavy quark energy loss. 
However, the discovery that the plasma produced in heavy ion collisions is itself a strongly 
coupled liquid has raised the question of how to understand the real-time dynamics of heavy 
quarks in a strongly coupled non-Abelian plasma. Treating all aspects of the dynamics as 
strongly coupled is of value hrst as a benchmark and second because it means that rigorous 
calculations of novel effects become tractable in plasmas with a gravitational dual. 

The simplest plasma in which one can calculate the rate of energy loss of a heavy quark 
is that in strongly coupled AA = 4 supersymmetric Yang-Mills (SYM) theory in the large 
number of colors (large Nc) limit, whose plasma with temperature T is dual to classical grav¬ 
ity in a 4-|-1-dimensional spacetime that contains a 3-|-l-dimensional horizon with Hawking 
temperature T and that is asymptotically Anti-de Sitter (AdS) spacetime [2, 3]. In the dual 
gravitational description, the heavy quark is represented by a string moving through the 
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AdS black hole spacetime, trailing behind its endpoint that follows the trajectory of the in¬ 
finitely heavy quark along the boundary of the AdS [4-7]. The earliest work on heavy quark 
dynamics in the equilibrium plasma of strongly coupled AA = 4 SYM theory [5-7] yielded 
determinations of the drag force felt by a heavy quark moving through the static plasma and 
the diffusion constant that governs the subsequent diffusion of the heavy quark once its initial 
motion relative to the static fluid has been lost due to drag. This work has been generalized 
in many directions since then. We will in particular need the modifications of the spacetime 
metric that describe a flowing, hydrodynamic, plasma in which there are gradients of the fluid 
properties as a function of space and time [8-10] . The corresponding modifications of the drag 
force were worked out to leading order in the fluid gradients in Ref. [11]. However, in this 
calculation the possibility of a nonzero density of some fermion species, and a corresponding 
nonzero chemical potential, was not taken into account. 

We begin in Section 2 by introducing the dual gravitational description of a strongly 
coupled plasma with both a chemical potential and fluid gradients, working to leading non¬ 
trivial order in both /u/T and the fluid gradients. At the same time, we introduce the dual 
gravitational description of the axial anomaly, relevant if the chemical potential is either that 
for left-handed quarks or for right-handed quarks. In Section 3 we turn off the anomaly, for 
example as appropriate if ^ is the chemical potential for baryon number, and calculate the 
corrections to the drag force in powers of /r/T, working to first order in fluid gradients. This 
yields the (straightforward although laborious) extension of the results of Ref. [11] to the case 
of a plasma with nonzero fi. The results of this Section constitute quantitative modifications 
to the drag force, but they do not introduce qualitatively new effects. 

In Section 4 we analyze a chiral plasma. That is, we take /x to be the chemical potential for 
either left-handed or right-handed quarks, and turn on the anomaly. In the dual gravitational 
theory this means we turn on the Chern-Simons term in the holographic action that we 
introduced in Section 2. This term gives rise to novel chiral contributions (contributions that 
change sign if the plasma contains right-handed quarks as opposed to left-handed quarks) to 
the charge and entropy currents and to the stress-energy tensor [9, 12], reproducing the chiral 
magnetic effect (CME) and chiral vortical effect (CVE) that had already been introduced 
without the use of holography [12-15]. These anomalous contributions to the hydrodynamic 
motion of a chiral fluid have been discussed widely [9, 12, 16-43]. These effects originate in 
topological aspects of the gauge theory [12, 14, 15, 33, 36, 37], as is of course the case for the 
axial anomaly itself [44]. 

All the previously analyzed consequences of the axial anomaly in a chiral plasma — the 
CME and its cousins — concern the motion of light quarks; in fact the quarks are usually 
assumed to be massless. We show in Section 4 that there are anomalous contributions to the 
drag force on an infinitely heavy quark that finds itself in a chiral plasma. The standard CME 
and related effects involve the generation of anomalous currents parallel or anti-parallel to 
an external magnetic field or the angular velocity vector of the fluid [12-15]. The anomalous 
contributions to the drag force that we compute, order by order in /x/T, have the same feature. 
This means that in the presence of a nonzero density of heavy quarks that are initially at rest 
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in the local fluid rest frame, the chiral drag force that we compute can yield a new contribution 
to the CME electric current, even though the heavy quarks themselves do not participate in 
the CME. Note, however, that if the plasma features equal and opposite number densities of 
some heavy quark and its antiquark, or densities of two species of heavy quarks with opposite 
electric charges, all the heavy quarks and antiquarks feel a chiral drag force acting in the same 
direction, and no electric current is generated. For example, the fluid produced in a heavy ion 
collision is seeded with equal numbers of charm and anti-charm quarks, and equal numbers 
of bottom and anti-bottom quarks. Since in any volume of the plasma in which there is an 
excess of, say, right-handed light quarks all the heavy quarks and antiquarks feel a chiral drag 
force in the same direction, the chiral drag force does not result in an electric current. Of 
course, if there were a nonzero chemical potential for some species of heavy quark, meaning 
an excess of those heavy quarks relative to their antiquarks, the push on all heavy quarks 
and antiquarks from the chiral drag force would result in a heavy quark contribution to the 
electric current. This would be an example of a correction to the CME or CVE currents, as 
has been found in other contexts [35, 38, 40-43]. 

It at first seems odd to find a nonvanishing chiral drag force on a heavy quark even when 
the heavy quark is at rest in the local fluid rest frame. We show in Section 5 that (as long as we 
neglect the gravitational anomaly) the resolution is related to the previously known fact that 
in the local fluid rest frame there are anomalous contributions to the entropy current, since 
if we go instead to a frame in which the local entropy current vanishes at the location of the 
heavy quark we find no chiral drag force on the heavy quark. In this frame, the heavy quark 
at rest is immersed in a flowing fluid, with nonzero momentum and charge currents that are 
both proportional to the axial anomaly coefficient, but the heavy quark feels no drag force. If 
we think of the heavy quark as a defect placed in these propagating anomaly-induced streams, 
the fact that the heavy quark feels no drag force is a direct consequence of the dissipationless 
nature of the CME current. The nondissipative character of the chiral magnetic effect has 
been discussed before [24]. Our calculation of the chiral drag force provides direct evidence 
for this fundamental attribute of the chiral magnetic effect. 

In Section 6 we discuss possible phenomenological consequences of the chiral drag force. 
The basic effect is that the heavy quarks and antiquarks in a heavy ion collision in which 
there has been a fluctuation resulting in an excess of right-handed (left-handed) light quarks 
will feel a force that pushes them in a direction perpendicular to the reaction plane that is 
antiparallel (parallel) to the direction of any magnetic field or fluid angular velocity vector. 
We shall show that the effects are small. Furthermore, they will average out in an ensemble 
of events. And, mesons containing heavy quarks are not so numerous in any single heavy ion 
collision. For all these reasons, it will be difficult to separate the effects of a small force that 
acts in the same direction on all the heavy quarks and antiquarks in one event from random 
forces that act differently on different heavy quarks. Perhaps clever correlation observables 
can be found. One possibility is to utilize the fact that in each event the direction of the 
kick that all the heavy quarks and antiquarks in that event receive, regardless of their electric 
charge, is opposite to the direction of the CME electric current in the light quark sector. 
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2 Holographic fluid 


To compute the drag force on a heavy quark moving through the strongly coupled plasma 
of AA = 4 SYM theory with a nonzero chemical potential and with spatial and temporal 
gradients in the fluid velocity, temperature and chemical potential we need to start with the 
perturbations to the dual gravitational theory that correspond to hydrodynamic flow [9, 10]. 
We shall present the corresponding 4+1-dimensional bulk metric in this Section, delaying the 
introduction of the heavy quark that we are interested in to Section 3. The dual gravitational 
theory is described by the 4+1-dimensional Einstein-Maxwell action 

* = - - 4^ / (2-1) 

whose equations of motion are 

Rmn + 4:Gmn + i^Fm^Rkn + = 0 

dN + ne^^^^^FNoFpQ = 0 . ( 2 . 2 ) 

Here, G 5 is the 4+1-dimensional Newton constant which according to the holographic dic¬ 
tionary, see e.g. Ref. [9], is related to the number of colors in the boundary gauge theory 
by Gs = Furthermore, Am and Fmn are a 4+1-dimensional vector potential and the 

corresponding field strength and k is the 4+1-dimensional Chern-Simons coupling, which is 
dual to the axial anomaly coefficient in the boundary gauge theory. In the case of A/" = 4 
SYM theory at strong coupling, 

(2.3) 

(See e.g. Refs. [12, 45, 46]). We should mention that we have chosen units in which the AdS 
radius RAdS bas been set to 1, meaning that our holographic dictionary coincides with that 
in Ref. [9]. The equations of motion (2.2) have a static black hole solution that corresponds, 
in the boundary theory, to a plasma with some nonzero temperature T and some nonzero 
density of right-handed fermions with chemical potential fi. Note that k is the bulk Chern- 
Simons coupling constant and it gives the boundary theory the axial anomaly. If we want a 
model for a theory like QCD in which there are both left- and right-handed quarks, we need 
to introduce two bulk vector potentials, one with k = — corresponding to jiR and the 
other with k = + corresponding to fip- If we then wish to consider only the case /tl = /i_R, 
we could define fiy = (mr + 1*l)/ 2, set the two bulk vector potentials equal to each other, and 
ignore k. Of course, k becomes relevant in any circumstance in which = (+i? — 7 ^ 0. 

In Eddington-Finkelstein coordinates, the metric and bulk gauge field that describe a 
static black hole of mass M and charge Q take the form 


1 


ds‘^ 


—r‘^f{r)u^Uydx^dx'^ + r'^ P^iudx^dx^ — 2u^dx^dr 
\/3Q 

o ^LL 


(2.4) 
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with 


(2.5) 


where is a constant vector satisfying = —1, where = r]^^, + and where 

we are working in axial gauge A^. = Q. The energy density, charge density, entropy density, 
temperature and chemical potential of the boundary theory strongly coupled plasma can be 
related to the M and Q of the dual black hole as follows [45, 46]. fi and T are given by 


= 


\/3Qi?AdS 


T = 


r+ 

^'^■^AdS 


2 - 


r^Y 

r+J 


r^y 

r+J 


( 2 . 6 ) 


where r_|_ and r_ are the larger and smaller real solutions of the equation /(r) = 0, and where 
we have temporarily restored the factors of i?AdS- These relations can be rewritten as 

/ \ 

:) 



p 2 

-^AdS 


37r2r2 


and 


2—^2 
r_ - 


-1 + 

1 


9 - 


16 


l + t/l + 


2^2 


j 


(2.7) 

The energy density, charge density and entropy density are given in terms of ^ and T by [9] 

-27rT , s = ^_+ 2.8 


e = 


3iV2 ^3 


, 

' B? 

■^AdS 




where Nc is the rank of the gauge group and we are working throughout in the large-A^c 
limit. The boundary theory plasma is conformal, meaning that it has Tji = 0, and so has an 
equation of state e = 3P, and has bulk viscosity C = 0- Its shear viscosity can be calculated 
via gauge/gravity duality and is given by j ^ [47, 48]. 

Analogously to the way that hydrodynamics is usually derived, the way to find a bulk 
metric that is the dual gravitational description of a flowing strongly coupled plasma is to 
look for a solution to the bulk Einstein-Maxwell equations in which T, /r and are all slowly 
varying functions of space and time, and to organize the calculation via a gradient expansion. 
The metric is expanded in powers of boundary gradients 


G;7. = G'W + GW + 0(a2), (2.9) 

where G® is defined in (2.4) and where gJJJ contains all possible gradient structures of first 
order, with unknown coefficient functions. The solution has been obtained up to first order 
in gradients and is given by [9] 

ds‘^ = —yf{r)u^j,Uydx^dx^ + yP^ydx^dx'^ — 2u^dx^dr + ^ F{r)a^i,dx^dx’^ 

+r‘^ja {P^Uy + PyU^) dx^dx’^ + ‘Jrr{d • u)Ufj,Uydx^dx^ (2-10) 

o 

where 

= P^^»py (9„n^ + dfiUo) - Ip^’^d • u (2.11) 

O 
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and 


1/ 2y/3Q^K„ T/ NO 

Ja — “I Mr^ 


( 2 . 12 ) 


where = efj,i,ai 3 u'^d°‘u^ is the four-vector containing the vorticity of the fluid, and where 
the functions F{r) and J(r) are, up to leading nontrivial order in fj,, given by 


F{r) = 


1 


47rr 

+ 


, 7rT\ 

2 arctan ( — I — log 


(r-b 7rr)2(r2 7r2r2) 




247r2r2 
6 


+ 


ttT 


3 47rT 

-^ + 

1 F 

arctan 


+ 


4(r -b ttT) 


VttT 


r -b ttT r2 -b 7r2T2 
^2 


-blog 


r2 -b 7r2r2 


+ O (/) 


(2.13) 


and 
J(r) = 




27r2rV(3r - 27rr) - 3(r^ - tt^T^) 


,4-Tn4\ 


247r3r4r2 

+0 (/r^) . 


, 7rT\ / (^ + ttT)^ \ 

2 arctan ( — ) — log ( ^ _‘)m 9 ) 


r J 


r2 -b 7r2r2 

(2.14) 


Note that since Q (x fj, the anomalous term in jo- that is proportional to Ki^ is of order 
while the leading ^-dependence in the metric (2.10) coming from F(r) and J(r) is of order 
fj?. We will see the effects of the terms of order in Section 3 and shall introduce the effects 
of order that originate from the axial anomaly in Section 4. In Section 4 we shall also 
introduce effects that originate from the axial anomaly that are proportional to an external 
magnetic field, rather than to the fluid vorticity. Doing so will require adding an additional 
term in j^-, a term that is of order and that is proportional to k and the magnetic field. 
Note that, as in Ref. [9], we are using the Landau definition of the fluid velocity u^, such 
that the Lorentz frame in which = (l,(l) corresponds to the frame in which the fluid 
momentum vanishes. In this local fluid rest frame, there can still be nonzero charge 
currents or entropy currents, as we shall see later. 

To the same order in /r/T and gradients as above, the bulk gauge field now takes the 
form 


.4 --^u 

■^(1 - O 


M + j^j.2 


(2.15) 


where we are still in = 0 gauge and where a(r) is a known function that can be found in 
Ref. [9]. However, the bulk gauge field does not affect the string that we will introduce in the 
next Section that constitutes the dual description of the heavy quark, and for this reason we 
do not require for our considerations. 

In most contexts relevant to heavy ion collisions, /x/T is small. If /x represents an excess 
of fiR over hl or vice versa, these only arise due to fluctuations or anomalous effects and are 
certainly small compared to T. If, as in Section 3, we treat a nonchiral plasma and think of 
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fj, as representing the chemical potential for quark number (which in turn is one-third that 
for baryon number), /r/T can be as large as ~ 1 in the lowest energy heavy ion collisions 
possible at RHIC but in higher energy collisions it is substantially smaller than 1. We will 
see in Section 6 that the effects of interest are suppressed by powers of /U^/(7rT) or /iy/(7rT); 
even if the quark number chemical potential pLy can get as large as ~ T it is always small 
compared to vrT. We will therefore work perturbatively in powers of pl/T. 

3 Drag Force 

In this Section we shall only consider p-i = fJ-R and shall therefore set k = 0 and think of pt 
as representing the quark number chemical potential {pii + 

The drag force has been calculated in a static plasma with pi = 0 and no gradients in 
Refs. [5-7]. The basic picture of heavy quark dynamics that emerges, with all but the initially 
most energetic heavy quarks being rapidly slowed by drag and then becoming tracers diffus¬ 
ing within the (moving) fluid, is qualitatively consistent with early experimental investiga¬ 
tions [49] . For a review, see Ref. [50] . Subsequently, the holographic calculational techniques 
were generalized to any static plasmas whose gravitational dual has a 4-|-1-dimensional metric 
that depends only on the holographic (i.e. ‘radial’) coordinate in Ref. [51] and heavy quark 
energy loss and diffusion has by now been investigated in the equilibrium plasmas of many 
gauge theories with gravitational duals [52-66]. In particular, the drag force on a heavy quark 
in a static plasma with pi 0 was calculated in Refs. [51, 52]. We reproduce these results 
in Section 3.1. More recently, the drag force on a heavy quark moving through the far-from- 
equilibrium matter present just after the collision of two sheets of energy density [67] has 
been calculated and compared to that in static strongly coupled plasma in equilibrium [68]. 
Motivated initially by the need to understand these results, in Ref. [11] the drag force was 
computed in a fluid with = 0 in which there are spatial and temporal gradients in the fluid 
temperature and flow velocity, to leading order in these fluid gradients. In Section 3.2, we 
shall reproduce these results and shall compute the leading order effects of fluid gradients in 
a strongly coupled AA = 4 SYM plasma with pi 0- (As an aside, we note that the effects of 
fluid gradients on photon emission have also been investigated [69].) 

3.1 Fluid at rest 

According to Refs. [5, 6], the dual picture of an (infinitely) heavy quark or antiquark moving 
through the strongly coupled plasma of W = 4 SYM theory corresponds to a string whose 
endpoint moves on the boundary of AdS along the trajectory of the heavy quark or antiquark 
of interest. The string trails “downward” into the bulk spacetime, toward the black hole 
horizon. It also trails behind the moving quark or antiquark. The drag force is defined as 
the force needed to move the heavy quark with constant velocity v through the plasma, so 
we take the endpoint of the string to move along the AdS boundary with constant velocity v. 
The drag force is determined in the dual gravitational description by the shape of the trailing 
string. Because the string trailing behind a heavy quark or a heavy antiquark moving with 
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the same velocity through the same fluid has the same shape, they feel the same force. This 
conclusion applies throughout: the introduction of gradients in the fluid, magnetic fields, and 
anomalous effects later in this Section and in Section 4 does not change the fact that heavy 
quarks and antiquarks feel the same force. We shall therefore henceforth refer only to heavy 
quarks. 

Without loss of generality one can consider the heavy quark velocity v to be directed 
along the x-direction and because we are considering a static plasma it makes sense to work 
in the rest frame of the medium, generalizing this only later. In this setup the probe string 
is described by the Nambu-Goto action 



where A is the ’t Hooft coupling, which we are assuming is large, and with g{T, a) the de¬ 
terminant of the induced worldsheet metric. The worldsheet metric is defined as gab = 
GAsdaX^dbX^ with x^{T,a) the string profile. Here, upper case Latin indices correspond to 
4-|-1-dimensional spacetime indices and lower case indices correspond to the 1-|-1-dimensional 
worldsheet metric indices. We also have freedom in the choice of the parametrization of the 
worldsheet and it is convenient to take 


t(r, a) = T , r{T, a) = a . 

For a quark whose trajectory in the boundary theory is x = vt, the ansatz for the shape of 
the string in the bulk that corresponds to the word “trailing” is 


x(r, a) =VT + ^(cr). 


(3.2) 


with y(r, a) = z{t, a) = 0. 

The string profile ^(fi) is found by extremizing the Nambu-Goto action and solving the 
resulting Euler-Lagrange equation 


dr 



+ da 




(3.3) 


We shall solve this for the string profile order-by-order in powers of g/T. The g = 0 (or 
equivalently zero quark number) solution was obtained in Refs. [5, 6]. We can write the 
expansion about the g = 0 solution in Eddington-Finkelstein coordinates as 

, vr — 2 arctan (^) 

27tT 


fO.2 = + o (/) 


(3.4) 


where the 0,2 superscript refers to zeroth order in gradients and second order in g/T and 
where the direction of is the same as the direction of v, which we are taking to be the x 
direction. From the Euler-Lagrange equation (3.3), we find that the equation of motion for 
^S’^(u) takes the form 


da{(J^ - 


7r^TS^)da^^’^{a) = RRS 


(3.5) 






where the left-hand side that we have given explicitly has the same form also at higher orders 
in /u/T while the right-hand side, an expression involving no derivatives with respect to a, is 
different at each order. Since the string endpoint is constrained to move along the constant 
velocity trajectory x = vt, we require CS’^(oo) = 0. Upon solving the equation of motion, we 
find that the string profile to order fj?/T'^ and to zeroth order in gradients takes the form 


X = VT + 


v{tt — 2 arctan (^)) 




27rr 




Utt^T \T 




4T(T 


5/2 


7 


-|-6 arctan 


7^17 


7r2r2 + (t2 


( ttt ) ' 72^/7 




arctan 


a 


ttT^J 


(3.6) 


where 7 = — v^. In the zero chemical potential limit, (3.6) is the solution obtained 

in Refs. [5, 6 ]. In solving the equations of motion to obtain (3.6), we fixed the integration 
constants by requiring that the worldsheet must be regular at the worldsheet horizon, to 
zeroth order in /t located at a = vrTya, and by requiring that the string endpoint follows the 
desired trajectory, as already mentioned. 

With the string profile in hand, we can now compute the drag force. The drag force is 
defined as (see e.g. Refs. [5, 6 , II]) 


F{t) = - hm v^^U;{T,a) 

rr—^ 


(3.7) 


where the string momentum is 


n- = [gradrX^ - Qrrd.X^] (3.8) 

and we are interested in the instantaneous drag force at r = 0. (Note that is not a Lorentz 
4-vector but one could introduce the proper force 7 /^ which has appropriate transforma¬ 
tion properties.) Upon substituting the string profile (3.6) into the definition of the string 
momentum we have 


lim n;)(o,cr) = - 


r2^ 

127r 


67 r 27 - (1 - 37 ) (3.9) 


from which we find that the drag force is given by 

r2^/A I 


- 

7(0) - 


127r 72 


{- iw ^ - < 5 ^) - (1 - 37) ^ 


(3.10) 


where = 7 ( 1 , u) is the four velocity of the quark. This result was first obtained in 
Ref. [51]. Note that this force, conventionally called the drag force, is defined as the force 
that an external agent must exert on the heavy quark in order to keep it moving at constant 
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velocity v. We see from the result (3.10) that / points in the same direction as v. The force 
that the fluid itself exerts on the heavy quark is —/, in the —v direction. 

We have obtained the drag force in the fluid rest frame, where = (1,0). We can now 
Lorentz-transform 7 /^ to a general frame, and upon doing so we find 



2tt 7 


{sw^ + u^) 



(1 + 3s) 
6 s 



(3.11) 


where s = u-w is the only Lorentz scalar that can arise at zeroth order in gradients. It also can 
be shown by direct calculation, starting from the definitions of the drag force and the string 
momentum, that w^f^{T) = 0. Note that in the nonrelativistic limit 7 —)• 1 we have s = —1 
and the ^-dependence of the drag force is ~ (l -|- whereas in the ultra-relativistic 

limit 7 00 we have s = —00 and the /i-dependence of the drag force is ~ (l -|- 

This means that in the presence of a nonzero chemical potential, the force required to move 
the heavy quark through the strongly coupled plasma is strictly speaking no longer a drag 
force, since the magnitude of the force is no longer proportional to the momentum of the 
heavy quark. This is also the case in the presence of nonzero fluid gradients, as shown in 
Ref. [11]. By convention, we will nevertheless continue to refer to the force needed to move 
the heavy quark through the fluid as a drag force. As a side remark here, but a side remark 
that we will need in Eq. (5.14) in Section 5, note from (2.7) that r_|_ ~ -kT (l -|- g(:^)^), 
indicating that in the nonrelativistic limit the drag force is proportional to r^. 


3.2 Corrections to the Drag Force Due to Fluid Gradients 

The drag force was calculated to first order in fluid gradients at zero pL in Ref. [11]. We can 
follow the same logic. We must first obtain the string profile to first order in fluid gradients. 
Since the string profile at zeroth order in gradients that we obtained in (3.6) above contains r 
only at linear order one can make the following ansatz, which satisfies the equation of motion: 

fhn ^ ^ ^ _ (3_^2) 


Here, the expression for the string profile to zeroth order in fluid gradients in a general frame 
can easily be obtained from our previous results and is given by 

f°’°(r, a) = VT-^ {u^v - u) (^arctan . (3.13) 

vanishes, and can also easily be obtained from our previous results. It can then 
be shown that choosing with the medium derivative defined as Dt = 

dt + v'^di, cancels the r-dependence in the equation of motion. Next, we expand about the 
zeroth-order-in-gradients solution, see the r-dependent terms cancel against contributions 
introduced by /i^’”((t), and find that satisfies an equation of the form 

= RHS (3.14) 
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where the right-hand side depends only on a, not on derivatives of a, and where the right- 
hand side is the same for i = y and i = z but is different for i = x. We shall not provide 
the expressions for the right-hand side here as they are lengthy. In solving the equations of 
motion (3.14), the boundary conditions are the same as at zeroth order and again we require 
regularity at the worldsheet horizon. The full solution for ^’*^(cr) can be found in Ref. [11]. 
We have obtained but again these expressions are very lengthy and we will not give 

them here, focussing instead on presenting results for the drag force. 

To first order in gradients and zeroth order in n, we reproduce the gradient corrections 
to the drag force on the heavy quark that were first obtained in Ref. [11] and that take the 
form 


hm 


VAttT 

27r 7 



sd^^s 




+C 2 {s)U^{d ■ u) — \/^(u • d)U^ 


(3.15) 


where s = u‘^Wa as before and where we have defined 


ci(s) = 7r/2 — arctan (\/^) — 'n'TF^{s) 

C2(s) = + (1 + s^)ci(s)) 

and introduced the projector = Ufj, + Here, the function F^{s) is obtained by setting 
// = 0 in the function F defined in (2.13). 

Turning now to the corrections that are leading nonzero order in both gradients and /r/T, 
following the logic set out above we expand the equations of motion for the string profile in 
a double expansion in gradients and n/T, solve them, and obtain the force from the string 
profile as we did at zeroth order. We find no contribution that is first order in /r/T. To second 
order in /r/T and first order in gradients, the correction to the drag force is given by 


rM — 

Rl, 2 ) “ 


v/A 
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— € 3 ( 5 ) (^U^{wd)s — — U^C 4 {s){d ■ u) — ^(y; . 5 ) log 

+f/'‘ ( W , a) ^ 
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(-S): 


- TTU/ i 2 I\ 6s^-I-(37r - 10)s^-I-4 - TT \ 

+ -—- 9 ) 5 + 5 ^ 03 ( 5 ) +-^- {u ■ d)U^ 

{-s)2 \ {-s)-2 J 

+ (sc3(s) - ——^ ) (u;-9)t/'^-b2c5(s)9'"log^ 

\ (-«)" / 


(3.16) 
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where we have defined 


C3(s) 

C4(s) 

C5(s) 


{-s)l 


3(-s)2 _ 


4 — TT + 6 (—s) 2 — 6 s^ + 4(— s)2ci(s) + 3(— s)2 log ( 1 
2(3^/^ + s(2 + 3s(—1 + \/^ + s 


—s ( 8ci(s) — 3s(l + s^) log ( 1 — 


-1 + s 


'—s 


2^s 

-1 + s 


(3.17) 


4 - 3s(a/^ + 2s) + 6\/^(s^ - 1) ^ci(s) + log ^ 


To a large extent, all of these contributions to the drag force can be considered corrections to 
the // = 0 results (3.15) from Ref. [11]; they introduce quantitative changes to the gradient 
corrections to the drag force, but do not change the story in qualitative ways beyond intro¬ 
ducing d^{^/T) terms that describe the contributions to the drag force on the heavy quark 
due to gradients in the chemical potential, and the corresponding charge currents. 


4 Anomalous Contributions 

In the previous Section, we have computed the complete correction to the drag force to first 
order in fluid gradients up to order in the absence of n, which is to say in the absence of 
any chiral anomaly, as for example for the case where and the ^ in the previous 

Section represents {fj-L + liJ?)/2. In this setting, the next order contributions to the drag force 
would either be of order or second order in gradients, and as far as we can see would not 
introduce any qualitatively new effects. 

We now allow /i/j and fiL to differ and hence we turn on k, the Chern-Simons coupling 
in the bulk gravitational theory that describes the axial anomaly in the boundary gauge 
theory. Upon doing so, we shall find contributions to the drag force on a heavy quark that 
are qualitatively new, arising in two ways. First, we find new effects that arise at order 
because the contribution to the metric (2.10) proportional to k is proportional to These 
effects are proportional to and hence are first order in gradients. And, 

they only arise when n represents a chiral chemical potential like /tl or meaning that 
the effects we shall consider in this Section should be smaller in magnitude than those in the 
previous Section. Second, we find new effects that arise at order and are proportional to 
K that arise only in the presence of a magnetic field. 

We shall also describe the contributions to the drag force that are introduced when we 
include the gravitational anomaly in the boundary theory. 

As is by now well studied [9, 12, 19, 20, 35], introducing the Chern-Simons coupling in 
the bulk corresponds to introducing anomalous contributions to the hydrodynamic equations 
for a chiral plasma (a plasma with //l 7 ^ l^Ft) in the boundary theory, contributions that 
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arise because of the axial anomaly in the gauge theory [44], The corresponding anomalous 
transport phenomena have been studied at weak coupling [13, 15, 23] as well as at strong 
coupling [9, 12, 19, 20, 23, 30, 35, 70]. These anomalous effects can be summarized by noting 
that in a chiral fluid in the presence of an external magnetic field B or of nonzero fluid angular 
momentum the axial anomaly causes both vector and axial currents to flow [12, 13, 20, 23]: 


Jv{x) 

Ja{x) 


^B + 

27r2 


27r2 


B + 


V 27r2 



Q , 


(4.1) 


where fiy = (mr + AiL)/2 and fiA = (/Ur — are the vector and axial chemical potentials. 

The first and second terms in the vector current have been named the chiral magnetic effect 
(CME) and the chiral vortical effect (CVE), respectively. The expressions (4.1) are valid for 
a t/(l) gauge theory; in Section 5 we will generalize them as appropriate for J\f = 4 SYM 
theory. We note also that in (4.1) we have left out terms that are higher order in chemical 
potentials (order fj? in the CME and order in the CVE) that we will need, and introduce, 
later in (5.4). The contribution to Ja that is proportional to arises due to the gravitational 
anomaly (see e.g. Refs. [23, 71, 72]). 

We shall introduce the four-vector B^ = Note that in the limit in which 

the fluid velocity is nonrelativistic, and take the form = (0, B) and ifj, = (0, 211) in 
the local fluid rest frame, with B and 11 being the (externally applied) magnetic field and the 
local angular velocity of the fluid. 


4.1 Chiral Vortical Drag Force 


In this Section, we shall compute the anomalous contributions to the drag force required to 
pull a heavy quark through a chiral plasma, in which currents like (4.1) are flowing. We begin 
by setting the external magnetic field B = 0, considering only the effects of vorticity in the 
chiral fluid. We derive and solve equations of the form (3.14) and find that at order the 
string profile contains no term proportional to r and confirm that a contribution at order 
can only come from the presence of the Chern-Simons term in the bulk metric. Direct 
calculation then yields a contribution to the string profile given by 


_ ^kT ///\3 Tr'^ry-A7r^r272(l-g7)o-2 ,^0-4 

37r27VT/ (T^(7r2T27-|-Cj2)(7r2r2 _|_ 0-2^ 


(4.2) 


and a vorticity-induced contribution to the drag force in a generic frame given by 


i ^ + {£■ w)w^^ 

U il,3 - ^ 


(4.3) 


This chiral contribution is directly proportional to the anomalous coefficient k, as expected. 
In the calculation above, /x represents /xr. We now repeat the calculation for /xr, flipping the 
sign of K, and find that when /xr and /xr are both nonzero the result is obtained by replacing 
/x^ in (4.3) by /x^ - /x| = 6^^A^J‘v + 


- 13 - 







4.2 Chiral Magnetic Drag Force 


Next, we are interested in the anomalous contribution to the drag force due to the presence of 
a magnetic field. (We will not attempt the analysis for the more generic case of background 
electric and magnetic fields.) First, we need to modify the gravitational metric in a way 
that corresponds to introducing an external magnetic field in the boundary gauge theory, see 
e.g. Ref. [39]. This is done by adding a term to the static metric (2.4) that corresponds to 
adding a new term to ja in (2.12) given by 


where 


7r2r2 wr; 

^ _ 27r2T2 tt^T* ^ ^ 

Csir) — K T 1-2(1 — 


CB{r)B^, 

(4.4) 

vr^rn / r2 \ 

(4.5) 

]^°^Vr2+7r2r2j ' 


The full consideration of hydrodynamics with external electric and magnetic fields can be 
found in Ref. [39]. We shall follow the same procedure as above. And, as we are interested 
only in the effects due to the presence of the magnetic field we can work to zeroth order in 
fluid gradients. Since there is no contribution in it that is linear in /r/T the string equation 
of motion (3.14) simplifies as in the vortical case above and we find a correction to the string 
prohle given by 



K //i\2cT2(7r2r272 + o-2)CB(CT)-7r2r272(7r2r2-hfj2)C'ij(7rr^)^ 

7r4r2 vr; (tt^TV - + ^ 2) ’ ^ ^ ^ 


where we are considering the magnetic field itself to be first order in gradients. This results 
in a magnetic-field-induced correction to the drag force given by 


(/^)i,2 = (^) s'^ CB{7^T^/^){B^ + {B■w)w^) . (4.7) 

We have excluded the term corresponding to the Lorentz force on the quark, as this does not 
depend on the medium or on the chiral anomaly and is not of interest to us here. 

We then repeat the calculation for flipping the sign of k, and hence replace fj? in 
(4.7) by ~ Bl = 4/i^/iy. Note also that we have assumed that the magnetic field is a 
perturbation; it could be interesting to perform a similar calculation in the case of a strong 
magnetic field. 

We defer speculating about possible phenomenological consequences of the vortical and 
magnetic contributions to the drag force that we have found in (4.3) and (4.7) to Section 6. 
We note, however, that both effects are small both because they are proportional to which 
arises only due to topological fluctuations in the plasma, and because they are proportional to 
fjLv/T (chiral magnetic drag force) or (chiral vortical drag force). This suppression is 

less severe in lower energy heavy ion collisions in which the quark number chemical potential 
is higher, but of course the lower the collision energy the rarer it is to produce heavy quarks. 
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4.3 Contributions to the Drag Force arising from the Gravitational Anomaly 

Before we continue, we introduce a further generalization. The term in the anomalous 
axial current in (4.1) raises the question of whether there is an analogue of this term in the 
drag force. Such a term could be substantially greater in magnitude in any phenomenological 
consideration, since is greater than in heavy ion collisions, and is much greater than 
The contribution to Ja in (4.1) is connected to the presence of a gravitational Chern- 
Simons term in the bulk metric (see e.g. Ref. [23]), namely a contribution to the bulk action 
that we have neglected up to this point that takes the form 


<55 = I d^xV^e^^PQ^AMRUpRAQR ■ ( 4 . 8 ) 


Here, Kg hxes the coefficient in front of the gravitational axial anomaly in the boundary 
theory. If, following Ref. [26], we add a single chiral fermion transforming under U{1)l to the 
holographic theory, we then have 



(4.9) 


We will use this expression when we make estimates. The value of the ratio (4.9) for the axial 
quark number current in two-flavor QCD is Kq/k = 3/20 while in the three-flavor case it is 
Kg/n = 3/16. Introducing Kg ^ ^ adds two new terms to the expression for in (2.12) that 
appears in the bulk metric (2.10), one in the direction of the fluid angular momentum and 
one in the direction of the external magnetic field. Working to zeroth order in gradients and 
in each case keeping the lowest two nonzero terms in the expansion in /r/T we find 


I 




4;u2 ^ 23 . 2)2 

7r2r2r8^ 


2 37r^r^\ (r^-F tt^T^) 

^ ) + vr2r2 (r2 - 7r2r2) \r^ + tt^T^ J 




3g,a = >^g^<T 


4^ - vr^r^T^ + 2r^ (r^ - vr^T^) log 


7r2r8T2 


Kg£cr 


( 3 (hr® - 


27r6 


7r6T6 


log 


r2 -K 7r2r2 


+ 0 (//^ 


1 


+ (I8r® - OvrV^T^ + vr^r^T^ - ISTr^r^T^ + Svr^T®) | + . (4.10) 


Note that Kg flips sign for left-handed quarks, which means that in a theory in which there 
are both left- and right-handed quarks the term in that is of order KgT^ will cancel, and 
we drop it henceforth. This means that the is only a quantitative correction to the in 
(4.4) that, given (4.9), is small in magnitude. In contrast, because includes a term that 
is linear in /i it can be substantially greater than the term in (2.12), as anticipated. 
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Completing the calculation, for the gravitational axial anomaly contribution to the drag 
force we find two terms, one in the B direction and one in the i direction, given by 

(4)f,i = (cb (r + {«■»)»") , (4.11) 

where we have only kept the lowest order terms in ^ that are nonzero in each case, and in 
particular have dropped the 0{^^) term in that could be obtained from the term 

in jg,a- 

The new contribution to the chiral magnetic drag force does not introduce any qualitative 
change. The new contribution to the chiral vortical drag force, however, is linear in /r/T 
whereas what we had found previously in (4.3) was of order meaning that the new 

contribution (4.11) is the leading one. At present no observational evidence of the term in 
(4.1) has been reported (except on the lattice, see Refs. [42, 73, 74])). So, if it is ever possible 
to see the chiral vortical drag force, in which (4.11) is the leading contribution, this could be 
a way to probe the physics of the gravitational axial anomaly in the gauge theory. 

5 The Dissipationless Character of the Chiral Magnetic and Vortical Ef¬ 
fects 

It was pointed out some time ago [24] using a symmetry argument that the charge currents 
induced by the chiral magnetic and chiral vortical effects must be dissipationless. The argu¬ 
ment in its essence is that these chiral effects are time-reversal invariant, and so cannot result 
in the production of entropy. In this Section we show that our results for the chiral drag force 
on a heavy quark provide an example of an explicit calculation that confirms the dissipation¬ 
less character of the CME and CVE. We do so by finding a particular setting in which we 
can place a heavy quark at rest within momentum and charge currents induced by either the 
CME or the CVE, much like placing a defect in a current-carrying wire or placing a rock in 
a flowing stream, and see that in the setting that we construct the CME- or CVE-induced 
current flows past the heavy quark defect without exerting any drag force upon it. 

Throughout this Section we shall largely neglect effects arising from the gravitational 
anomaly, in effect setting Kg = 0. In order to achieve our goals in this Section, we need 
not include spatial or temporal gradients in the fluid other than vorticity. Gradients can 
also be added, again at the expense of adding complication, and again without affecting the 
conclusion. 

We shall make our argument in three steps. First, we begin by noting that a heavy quark 
that is at rest in the local fluid rest frame feels a chiral drag force in a chiral fluid, in the 
presence of either a magnetic held B oi a huid vorticity 12. The local fluid rest frame is the 
frame in which the huid momentum vanishes at the location of the heavy quark. If the 
heavy quark is at rest in this frame, then in this frame = w°‘ = ( 1 , ()), 7 = 1 and s = — 1 
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and the leading contributions to the chiral drag force due to B and 0 that we have calculated 
in (4.7) and (4.3) simplify to 


Best- 27r3 r2 SttS r2 ’ ^ ’ 

where we have used the fact that Cb{t^T) = 1. Recall that jjL here represents When 
and Hr are both nonzero, is replaced by AhaHv and is replaced by QhaHv + 
striking that a heavy quark at rest in a stationary chiral fluid with ha > 0 {ha < 0) feels a 
force exerted on it by the chiral fluid in a direction antiparallel to (parallel to) the direction 
of a magnetic field or fluid vorticity. (Recall that k = — is negative for h = and recall 

that / is the force that some external agent must exert to keep the heavy quark at rest while 
—/is the force that the heavy quark feels from the fluid.) 

Second, we ask what we expect will happen if we release the heavy quark at rest and 
allow it to move under the action of the force —/, with / as in (5.1). The heavy quark will 
start to move through the fluid with some initially increasing velocity v. As long as v is small 
in magnitude we can neglect the resulting change in the chiral drag force (5.1), but because 
the heavy quark is now moving it will feel an ordinary drag force (3.10) (as calculated in 
Refs. [5, 6]) in addition, which for small v is given by 


/drag = • (5-2) 

This means that the heavy quark will accelerate until it reaches a terminal velocity 

h‘^B ^ 4k H^n 

^terminal - K ^ (Trr)"^ 

at which it experiences no further net force. We therefore reach the remarkable conclusion 
that in the local fluid rest frame a heavy quark at rest feels a force (5.1) while a heavy quark 
moving through the fluid with velocity ^terminal feels no force. ^ 

Third, we boost by a velocity ^terminal to a new frame in which the heavy quark that was 
moving with velocity ^terminal in the original frame is now at rest. In this new frame, the fluid 
at the location of the heavy quark is flowing with velocity —^terminal) the heavy quark is at 
rest, and the heavy quark feels no force. The fluid momentum current is proportional to k 
and is entirely due to the CME and CVE, and the heavy quark placed in the flowing fluid 
feels no force. This is an explicit demonstration of the dissipationless character of the chiral 
magnetic and vortical effects. 

In the remainder of this Section, we will generalize this conclusion in two steps. First, we 
will show that the frame in which a heavy quark at rest in the flowing fluid feels no force is, 

^Note that if there were a chemical potential for heavy quarks, meaning an excess density of heavy quarks 
compared to the density of heavy antiqnarks, then once all the heavy quarks and antiquarks are moving at 
■yterminai there would be a resulting heavy quark current in the local fluid rest frame. 
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in fact, the frame in which the local entropy current vanishes. We will then use holography 
to demonstrate that this conclusion generalizes. 

The full expressions for the charge current and the entropy current including the chiral 
magnetic and chiral vortical effects are [12, 13, 20, 23] 



^ { 1 nu' 

nu^ + C\^i- - — 

1 Ci?s 1 

- V- -Br - 


1 


B, + -C 



3e + Py 


(5.4) 

(5.5) 


where n is the number density of right- or left-handed quarks and s is the ordinary entropy 
density. In these expressions, C is the coefficient of the chiral anomaly, 9^ = CE ■ B, and 

is given by C = l/( 47 r^) in a 17(1) gauge theory as in (4.1) while here, in AA = 4 SYM theory, 
we have 


r = -^ = JfL_ 

47r2V3 ' 


(5.6) 


We see that when the fluid is at rest there are nonvanishing charge and entropy currents. 
However, if we boost to a new frame in which the new (primed) fluid velocity is related to 
the original by 


_ B - 

^ 2e + P ^ 2 ,e + P ^ 


, 1 c^^ 

= Ua - - 


(5.7) 


then in this frame the charge and entropy currents are given by 


+ C^iB^ + , (5.8) 

s'^ = su'^ , (5.9) 

meaning that we have boosted to a frame in which u'^ = ( 1 , 0 ) means no entropy current.^ 
Note, of course, that when u'^ = (1,0) there is a nonzero momentum current, 7 ^ 0; the fluid 
momentum vanishes when = (1,0). In the nonrelativistic limit, the boost velocity needed 
to accomplish (5.7), namely to boost from the local fluid rest frame to the local entropy rest 
frame, is 

Vhoost = 

where we have used the nonrelativistic approximation I ~ 2H. Upon using (5.6), noting that 
£ + P = /2, and comparing to (5.3) we see that 


^^boost — 'i^terminal • (^-H) 

We conclude that the frame in which a heavy quark at rest in the flowing fluid feels no force 
is in fact the local entropy rest frame. 

^In the presence of the gravitational anomaly, i.e. if Kg 7^ 0, the entropy current (5.5) includes additional 
terms and, consequently, does not vanish in the “no drag frame” [75]. 
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The generality of this conclusion will become apparent after it is recast holographically, 
in terms of the dual gravitational description of the fluid. In particular, upon so doing we 
will see that although in the derivation above we only used the contributions to the chiral 
magnetic and chiral vortical drag force and to the ordinary drag force that arise at the lowest 
nontrivial order in ^/T in each case, the conclusion that we have reached in fact holds to any 
order in /r/T. 

The authors of Refs. [28, 32, 76, 77] have shown that in the holographic description of a 
flowing strongly coupled plasma with = 0 with the 4+1-dimensional metric (2.10) the local 
entropy rest frame is the frame in which the metric function in (2.10) vanishes at r = rh, 
where Vh is the location of the outer horizon, namely the largest solution of /(r) = 0, with 
/(r) also being a metric function in (2.10). For the case of a static fluid with fi = 0, rh is given 
simply by rh = ttT. In a static fluid with /r 7^ 0, we have instead rh = r^ where r_|_ is given in 
terms of + and T in (2.7). We must therefore compute the force on a heavy quark at rest in 
the frame in which j{rh) = 0. This computation turns out to be sufficiently tractable that we 
can push it far enough without expanding in fi/T. We shall include a nonzero magnetic field 
and vorticity. As throughout this Section, we neglect the gravitational anomaly and fluid 
gradients. In this case, ja takes the form 

ja = W^ir)B^ + W^ir)i^ (5.12) 

where the effects of the chiral anomaly enter through the functions W^{r) and W^{r) whose 
form we obtained explicitly to leading nontrivial order in +/r in (4.4) and (2.12), respectively. 
The force on the heavy quark at rest can be calculated following the procedures developed in 
previous Sections, and one finds 

/ = (rlW^{rh)B + rlW^{rh)i). (5.13) 

It can be also verified that if one expands this answer for the anomalous drag force in powers 
of n/T, it does indeed coincide with (4.7) and (4.3). If we balance this force against the drag 
force on a slowly moving quark, which in this general context is given by 

/ = , (5.14) 

we see that the terminal velocity of the heavy quark is given by 

Winai = {w^{rh)B + W^{rh)i) . (5.15) 

Comparing this result with (5.12), we see that in the local entropy rest frame in which 
jai^h) = 0, the terminal velocity vanishes: ^terminal = 0. This means that in the local entropy 
rest frame, a heavy quark at rest feels no force. Note that in obtaining this result via this 
holographic calculation we did not need to expand W^{r) or W^{r) in powers of n/T. This 
means that the result (5.15) is valid to all orders in /r/T. Of course, in (4.4) and (2.12) we 
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have explicit expressions for VK^(r/j) and W^{rH) only to lowest nontrivial order in at 
higher orders, and W^{rH) receive corrections but the form of the expression (5.15) 

for the terminal velocity of the heavy quark remains unchanged. 

Note also that although we have not given the resulting more complicated expressions 
here, if one adds the effects of fluid gradients into the expressions for ja, the entropy current, 
and the drag force, while keeping Kg set to zero, the conclusion that a heavy quark at rest in 
the local entropy rest frame feels no force remains unchanged. 

Let us close this section by restating this general result in its simplest form. In a strongly 
coupled fluid in which ///j > in the presence of a magnetic field B and for simplicity in 
the absence of any fluid gradients including vorticity, we have analyzed the chiral drag force 
on heavy quarks with two different velocities: 

• If the first heavy quark is at rest in the local fluid rest frame, meaning that it sees 
around it a fluid with no momentum flow, this heavy quark feels the fluid around it 
exerting a chiral drag force on it pushing it toward the direction oi —B. It also sees 
around it a charge current in the direction of B and an entropy current in the opposite, 
—B, direction. (All signs are reversed if instead the fluid has > fiR.) 

• If a second heavy quark is released and allowed to accelerate under the influence of the 
chiral drag force, it does so until it reaches a terminal velocity in the —B direction. If 
we now boost to the rest frame of this heavy quark, we find a heavy quark that feels 
no force. This heavy quark sees a charge current in the B direction that is larger in 
magnitude than that seen by the first heavy quark. This heavy quark also sees the 
fluid around it flowing, with fluid momentum in the B direction. (Again, all signs are 
reversed if instead the fluid has fiL > hR-) 

The second heavy quark provides an example of a “defect” at rest in a flowing chiral fluid, 
with both fluid momentum and charge flowing past it, and yet feels no force. Thus, it provides 
an explicit example of the dissipationless character of the chiral magnetic and chiral vortical 
effects. 

This becomes particularly important now that the chiral magnetic effect has been seen [78] 
in a condensed matter system, namely the Dirac semi-metal ZrTes, since we can assume that 
the zirconium pentatelluride crystals used in the experiment contain some defects. 

6 Outlook and Potential Phenomenological Consequences 

In Section 3.2, we obtained the contributions to the drag force on a heavy quark that are 
second order in ///T and first order in fluid gradients in a fluid in which m = /r/j, meaning 
that there are no contributions to the drag force coming from the chiral anomaly. These 
results generalize the /i = 0 results of Ref. [11] to nonzero //. In Ref. [11], the analytic results 
for the contributions to the drag force that are first order in gradients were used to obtain 
an understanding of a variety of curious features of the drag force on a heavy quark caught 
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between two colliding sheets of energy density, features that were discovered numerically in 
the holographic calculation of Ref. [68]. We look forward to seeing our results from Section 3.2 
used similarly, in conjunction with some future calculation of the collision of two sheets of 
energy and charge density. 

We have calculated the anomalous contributions to the drag force on a heavy quark in 
a chiral plasma to lowest nontrivial order in jx/T in Section 4. These contributions to the 
drag force change sign in a chiral plasma with hr > jXL relative to that with fiL > HR- 
This parity-odd symmetry makes the chiral drag force novel, qualitatively different from all 
contributions to the drag force on heavy quarks that have been calculated previously. Its 
unique features — namely that it is a force that pushes all heavy quarks and antiquarks in 
the same direction, either parallel to or antiparallel to the magnetic field vector B or the 
fluid vorticity vector Q, with the direction of the force being opposite in a chiral plasma with 
Hr > Hl relative to that with hl > HR — make it detectable at least in principle via suitable 
observables defined for this purpose. In practice the effect is likely to be small in quantitative 
terms. We shall provide a rough estimate of the magnitude of the chiral drag force, before 
describing a possible observable. 

The full leading order expressions for the chiral drag force are given in Section 4 for a 
quark moving through a chiral plasma with any velocity. Here we recapitulate our result for 
the anomalous force on a heavy quark at rest in the local fluid rest frame 
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( 6 . 1 ) 


where we have restored the effects due to the gravitational anomaly, proportional to Hg, that 
we dropped in (5.1). We can compare this result for the chiral drag force on a heavy quark 
with u = 0 to the result for a heavy quark in the u —)• 1 limit, where we find 


27r 

71 


/ 


4 (k — 2Kg)H^ 

37r2r2 


{B 


32KgH,^ 
v)v - 


v)v — 


Ah^{3k — AOKg) 
97r2r2 


(H • v)v 


( 6 . 2 ) 


There are certainly differences between (6.2) and (6.1), for example the fact that in the 
ultrarelativistic limit the force is in the v direction rather than in the H or H direction, 
although its magnitude is greatest when v is parallel or antiparallel to H or H. The principal 
message coming from (6.2), however, is that the anomalous drag force on a heavy quark with 
u —?• 1 is comparable in magnitude to that on a heavy quark at rest. All heavy quarks and 
antiquarks feel a chiral drag force whose magnitude is only weakly dependent on their velocity. 
Given that, it seems clear that the fractional effects of the chiral drag force are largest when 
u ~ 0 rather than when u —)• 1. Although in principle there are anomalous forces on acting on 
the ultrarelativistic b-quarks that become b-jets, it seems a better bet to focus on b-quarks 
that are almost at rest relative to the fluid. We shall see that even in this case the effects are 
small. 

Let us consider a heavy quark with mass m that is initially at rest in the local fluid rest 
frame in a chiral plasma in the presence of a magnetic field B. This quark feels a chiral drag 
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As we described in Section 5, this heavy quark starts to move due to this force and as it moves 
it begins to feel the standard drag force (5.2) also, meaning that the heavy quark accelerates 
until it is moving with a terminal velocity 


^terminal /til 


16k. 




3 k J (vrT)' 


,B. 


(6.4) 


How long does it take for the heavy quark to accelerate from rest to a velocity that is close 
to ^terminal? This timescale is parametrically of order m|uterminai|/|/^| 27rm/ {y/X{7rT)^), 
meaning that for charm (bottom) quarks it is of order 1 fm/c (a few fm/c). For our rough 
purposes we can therefore estimate that the chiral drag force gives all heavy quarks a momen¬ 
tum that is around m ^terminal- Bottom quarks produced at rest in the fluid therefore pick up 
a contribution to their momenta arising from the chiral anomaly that is of order 


^terminal — ^terminal 





n'^B 


= —0.449 rufe 


fJ^VfLAB 

(7rr)4 


~ —3 MeV 


mb 


HV 
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(6.5) 


where we have used (2.3) and (4.9), and replaced by For concreteness, we have 

used the AA = 4 SYM values in evaluating the purely numerical factor k (1 -3^); QGD 

all components of this coefficient — namely k, Kq/k and the 16/3 — will take on different 
values. There will certainly also be other places where our AA = 4 SYM calculation differs 
from QGD by factors of order unity. Note, finally, that we have throughout used units in 
which a factor of e has been absorbed into B and E. Restoring this means replacing B by 
eB in (6.5). 

Reading the final expression (6.5) from left to right, we see that: 


• The effect is small. 


• The effect can be enhanced by lowering the collision energy, as doing so increases fiy, 
but one should not lower the collision energy too far since doing so makes heavy quark 
production rarer and also reduces the initial temperatures reached in the collision. Gen¬ 
erating the chiral effects we are discussing requires the presence of quark-gluon plasma. 
Recent measurements [79] indicate that the observables that are sensitive to the chiral 
magnetic effect that were seen previously in high energy heavy ion collisions [80-83] 
persist robustly down to collision energies corresponding to baryon chemical potentials 
/Us > 0.3 GeV, meaning /Uy >0.1 GeV. 
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• It is difficult to estimate the magnitude of the fluctuations in ha in heavy ion collisions. 
The authors of Ref. [84] have recently estimated that ha can locally be as large as 0.1 
GeV in heavy ion collisions at top RHIC energies. Estimating this more reliably will 
require the development of a relativistic viscous chiral magnetohydrodynamic code for 
heavy ion collisions. 

• In heavy ion collisions with a nonzero impact parameter, a magnetic field is created 
initially by the charged spectators. The conductivity of the quark-gluon plasma slows 
the decay of the magnetic field, delaying the decay of the magnetic field in the plasma 
to long after the spectators are gone. Reliable estimates here also await a future mag¬ 
netohydrodynamic analysis of heavy ion collisions, but perhaps \B\ can remain as large 
as (0.1 GeV)^ for a few fm/c [85-88], although this is more likely to be an overestimate 
than an underestimate. 

• Taking the average temperature seen by the heavy quark as vrE = 0.5 GeV is on the 
low side for a benchmark value with which to make an estimate, although maybe not 
unreasonable for heavy ion collisions with hv >0.1 GeV. 

Our estimates suggest that the chiral magnetic drag force on bottom (charm) quarks and 
antiquarks pushes them all to the extent that it gives them all a common momentum that 
might be as large as ~ 3 MeV (~ 1 MeV), although it should be clear that this estimate is 
first of all at present very crude and second of all is more likely to be an overestimate than an 
underestimate.^ Furthermore, in events in which there are some regions of the plasma with 
EJ? > El s-iid other regions with hr < ELj some heavy quarks will be pushed in the direction 
antiparallel to B while others will be pushed parallel to B, reducing the net effect in the event 
as a whole. 

Although the effects of the chiral drag force are small, we do not wish to underestimate 
the ingenuity of our experimentalist colleagues. Perhaps they can devise sufficiently sensitive 
correlation observables to see the small effects of the chiral drag force on heavy quarks. For 
example, it is worth constructing event-by-event observables that can see whether the B and 
D mesons in an event have picked up a net momentum perpendicular to the reaction plane in 
the “downward” (“upward”) direction in those events in which the chiral magnetic effect has 
resulted in an electric current “upwards” (“downwards”) with positive light hadrons pushed 
“upwards” (“downwards”) and negative light hadrons pushed “downwards” (“upwards”). A 
correlation observable like this uses the GME current to define the direction in which the 
heavy quarks should be pushed by the chiral drag force, and then checks whether a net push 
on all the heavy quarks and antiquarks in the event in this direction is seen. A nice feature 

^From (6.1) we see that there is a chance that the chiral vortical drag force could be larger than the 
chiral magnetic drag force since although Kg k the chiral vortical drag force is larger by some purely 
numerical factors and because it is not suppressed by hv/{t^T). However, at present the magnitude of H in 
the droplets of fluid produced in heavy ion collisions is poorly constrained and we have also only just seen the 
first measurements of observables that can receive a contribution from the chiral vortical effect [89]. We will 
therefore not attempt an estimate of the magnitude of the chiral vortical drag force. 
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of such a correlation observable is that the observable effects of both the CME current and 
the chiral drag force should each be equally suppressed by the partial cancellation between 
regions of the plasma with ^r> and regions with < /tl- 

Given the estimates that we have made, it is certainly our impression that the chiral drag 
force is of interest principally from a theoretical perspective, rather than as a phenomenon 
that experimentalists should expect to observe in heavy ion collisions. Even if they are only 
observable in principle, though, it is remarkable to see effects due to the chiral anomaly giving 
all the heavy quarks and antiquarks in a chiral plasma a kick. In Section 5 we have also used 
our calculation of the chiral drag force to provide explicit evidence for the dissipationless 
character of the chiral magnetic and chiral vortical effects, showing that the currents that 
they describe can ffow around defects without any dissipation. 
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